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We consider the micellization of block copolymers in solution, employing self consistent field 
theory with an additional constraint that permits the examination of intermediate structures. From 
the information for an isolated micelle (structure, binding energy, free energy) we describe how the 
global thermodynamics of these systems can be obtained, which can be used to build a realistic 
phase diagram; the role of translational entropy must be addressed in this regard. 



I. INTRODUCTION 



Amphiphilic molecules, presenting simultaneous solvophobic and solvophilic interactions, form many interesting 
structures in solution in which the solvophobic residues are screened from the solvent, such as bilayer structures 
and micelles. The later aggregate is remarkable, since, as they begin to proliferate around what is known as the 
critical micelle concentration (CMC), its size does not become macroscopic, but is limited by the particular molecular 
features. It is the main goal of this paper to provide a unifying view of micelle formation for the case of copolymers 
in solution!, in the regime where the micelle concentration is low, so that each micelle can be regarded as isolated; 
the translational entropy of the aggregates can play an important role in this regime, as we will discuss. (A dense 
fluid will of course be more complicated, since the interactions between micelles must be taken into account 2 ^.) 

We will use a simple version of self consistent field theory (SCFT) 4 ' 5 . The use of self-consistent theories for 
copolymer micelles has a long history, beginning with the work of Yuan et al. in 1992 6 . Our approach differs from 
previous ones in that we are not primarily interested in the calculation of properties of isolated micelles, but rather 
in how to include this information in a thermodynamic description of the micellar fluid. Even though the importance 
of the excess free energy of the micelles in defining a proper CMC and distinguishing between competing structures 
has been recognized since the first works!, we feel the connection between it and the concentration of micelles, which 
can lead to differences betweenpossible definitions of the CMC, has not been explored in detail. Additionally, we also 
employ the method of Refs.^lijin order to explore intermediate structures, that lie between equilibrium ones. 

This method can be classified with previous approaches to obtain the general thermodynamics from the properties 
of a single aggregate. There are several studies that focus on planar and curved bilayer— or monolayer— mem- 
branes, obtaining quantities like surface tensions and bending rigidities; micelles can be considered as limiting cases 
of these structures, as is discussed in some case a 11 ' 12 . There are studie s focused mainly on micellar structure and 
thermodynamic a 13 ' 14 , but few consider a global treatment. (Refs. ITHT^ come quite close, but are mainly concerned 
with periodic dense phases.) 

In section [n] we calculate properties of individual micelles; we discuss how to build the global thermodynamics of 
the system in section ITTT1 



II. ISOLATED MICELLES 



We have decided to focus on a copolymer-homopolymer mixture 1 ^, using the standard Gaussian model for polymers 
treated within SCFT, perhaps the simplest molecular theory that produces realistic results for micellization (just as it 
is, arguably, the simplest one that produces realistic complex mesophases). This way, we can avoid the complications 
associated with similar models that are more refined and specifioiSiii, while keeping a molecular theory that is closer 
to experimental systems than other models based on hard sphere fluidsiSiiiiiSii^. 

We consider an incompressible mixture of copolymers and homopolymer with volume V. The copolymers have N 
statistical units of volume p^ 1 and statistical length a; these are of type A from one end to monomer fN and type 
B from there on. The homopolymer is made of monomer A and has aN units of the same volume and length. In 
a system of volume V there will be an overall homopolymer volume fraction </>, which will be the spatial average of 
a homopolymer distribution 4>(r). By incompressibility, the corresponding copolymer volume fraction will be 1 — </>, 
and its distribution, 1 — 4>{r). Notice the overall copolymer concentration is p c = po(l — (f>)/N , and the homopolymer 
one, p h = p (f)/(aN). 

We will consider a grand free energy 

kTpoV = -( z <lh + q c ) + xN J dr<p A (r)<p B (r) + J dr [(f> A (r)w A (r) + (f>B(r)w B (r)] (1) 
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FIG. 1: Excess grand free energy as a function of r\ for values of (f> b — 0.61,0.65,0.70,0.75. The dotted lines span the range 
for which no solutions are found. The points correspond to the selected values that are used in Figure 2. 



The first term in the right side is the configurational entropy, and q c being the partition functions of a homopolymer 
and copolymer in the corresponding fields and z the fugacity of the homopolymer (by incompressibility, that of the 
copolymers can be taken to be 1, shifting the scale of chemical potentials). The second term is the interaction 
energy; the third contains the coupling with the fields wa and the fourth, a Lagrange multiplier, £(r), to enforce 
incompressibility. The last term introduces an additional variable ip which lets us sweep a range of metastable 
structures by "pinning" the profiles to a certain value at some point in space r = tq. (This goes beyond standard 
SCFT, as was used for this model in Ref. ^| even though this technique was used by the same author in Ref. for 
a similar system; see also Ref. for an application to related structures, spherical nucleation bubbles, and Ref. ITsl 
in which such a field is helpful to study bilayer fusion; in Ref. 19 hards wall are used in order to explore different 
micellar sizes.) 

The extremization of the grand free energy with respect to all the volume fraction profiles and fields leads to 
the self-consistent equations that have to be solved, with the copolymer and homopolymer profiles related to the 
fields through Green functions, g c (r, s), g^(r, s), which satisfy the standard diffusion equations^ and also provide 
the partition functions Vq c — J drq c (l), Vqt = J drqh(a). This kind of theory is well know to produce rich phase 
diagrams, with many periodic structures. But it is also possible to obtain certain structures that are localized. 
In these, the volume fraction profiles of the different components tend to their bulk values away from the spatial 
point where the structure is located. The bulk values are those of the corresponding disordered phase: <f)(r) — ► <f> b , 
4>a(t) — > 4> h + /(l — fc ), and <^s(r) — > (1 — /)(1 — <p b ). These localized structures correspond to isolated bilayers and 
micelles. 

The theory is simple to implement in planar, cylindrical and spherical geometry, but we will focus on the later 
in this article. The diffusion equations have been solved in real space using Cranck-Nicholson's method, which is 
much simpler conceptually (even if not so efficient computationally) than Fourier methods. (This procedure is also 
useful for problems in which the geometry is not knownS 1 or which are not very symmetric^.) With the help of the 
additional tp field, we study spherical micelles, as well as larger, metastable structures, corresponding to spherical 
bilayer vesicles. Notice this choice of tq as a "reaction coordinate" has some problems. Bilayer structures have two 
points at which 4>a{t) — 4>b {f)\ for spherical and cylindrical geometries, it is preferable to assign ro to the outer one 
(the one for which r is larger), since the inner one will be seen to disappear as the profiles become micellar. More 
importantly, for some values of ro, typically at regions of transition between different morphologies, no solution to 
the SCFT is found, a fact we will discuss below. 

We employ the iterative method described in Ref. I20I modified to include the ip field, to solve the self-consistent 
equations. It is found that 200 space points per ay/W and 4000 points along the chain are enough to provide results 
that cannot be distinguished from previously published results. A system size of 5ayN is sufficient to make finite 
size effects negligible. 

We will choose this set of parameters: %iV = 12, / = 0.6, a = 1; in Ref. we see this system shows, at high 
homopolymer concentrations, an "unbinding" of a micellar bec crystal into isolated spherical micelles. In Fig. 1 we 
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plot Afi, the excess grand free energy, as a function of the Tq for several values of (f> b , for spherical symmetry. We 
can see that for larger values of (j> b (i.e., low values of bulk copolymer), there are no nontrivial stable or metastable 
solutions. Above a certain value, there appears a minimum, which corresponds to a spherical micelle (this is the 
value associated with the CMC in some works). The Af2 is still positive, but becomes negative at lower <pb, a value 
that would be close to the CMC. Then, at even lower values, the asymptotic straight lines change slope. This linear 
behavior, for large aggregates, is characteristic of spherical vesicles, and, as described below, the corresponding profiles 
are indeed vesicular. The free energy as a function of the area A will tend to a lin o 8 ! 11 : AQ — > 7 A + 8it(k + Kg), with 
slope 7, the surface tension; n is the bending rigidity and kq is the Gaussian bending rigidity. That is, the change 
in slope corresponds to a change in the surface tension of the membrane; indeed, the condition for the stability of 
a membrane is 7 = 0. There usually is some leeway regarding the definition of the area A (i.e., we could propose 
A = 47rrQ, hence our choice of x axis), but when 7 = there is none. This method can be used to obtain the rigidities, 
since the cylindrical geometry can separately provide k 810 . In this particular case, though, the lamellar transition is 
preempted, since at this point the AO of the micelles is quite negative, as is clearly seen in the curve for <j>h = 0.65 in 
Fig. 1. This is just what would be expected, since our CMC is just the "cubic phase unbinding" of Ref. H, i.e., the 
stable periodic phase corresponding to our choice of parameters is a cubic phase, not a lamellar one. 

In Fig. 2 we present some typical profiles for the case cf> b — 0.65, at the points in Fig. 1. Fig. 2(a) corresponds to 
the stable micellar structure, with a core composed of B copolymer and a corona of A copolymer. As tq is increased 
the micelle grows, Fig. 2(b) shows the most "swollen" micelle we are able to obtain, still with the same architecture. 
After a region of ro for which no stable profiles are found, the next stable structures show a core composed of the A 
region of the copolymer and a little solvent, Fig. 2(c). This sort of "proto- vesicle" is a micelle with A-B-A structure, 
instead of the previous B-A one. For larger values of ro homopolymer progressively fills the core (which is now of an 
A nature, and thus compatible with it), while the copolymer profiles tend to the planar A-B-A bilayer structure as in 
Fig. 2(d), which corresponds to a local maximum in the excess free energy. 

We would like to point out it is possible to connect this results with an approximation proposed in Ref. Il2l to 
obtain the whole set of curves for unstable structures from information about critical ones (i.e., at local minima or 
maxima). Our results show that this approximation works quite well for large radii, but breaks down at small ones. 
In fact, an examination of the structures shows that the collapse of micellar structures at lower values of <fi b does not 
take place because they become unstable with respect to planar membranes (as is the case in Ref. Il2|) . but because 
the whole system begins to develop long range oscillations — that is, because we reach the spinodal. 
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III. GLOBAL THERMODYNAMICS 

We have seen one of the most important micellar quantities is its excess (grand) free energy, AJ7; another one is the 
excess solvent, A<p, a negative quantity (the excess of copolymer is — A<f>, positive). At low micelle concentrations, the 
components are either in the bulk or in micelles, so the total concentration of copolymer p t c — p h c + p m AN c , where p h c 
is the bulk concentration, p\ — po(l — (f) b )/N, AN C is the excess of copolymer molecules in the micelle (proportional 
to \A</>\, as we discuss below) and p m is the concentration of micelles. This quantity is expected to satisfy^ 

p m = — exp{-ASl/kT), (3) 
v m 

which is dominated by the condition AQ = 0, but the value of the volume v m , which is related to the translational 
entropy of the micelles, can make a difference, leading to different choices for the CMC. (Another important discrep- 
ancy is that many authors choose to relate the CMC with the point at which the micelles become stable (in general, 
metastable) 12 i 14 i 16 i 21 ' 22 .) We propose, based on our previous worki^: v m = (po/N) exp[(AE — AF)/(kTAN)], where 
AE is the excess interaction energy, (second term in Eq. ^ minus the bulk energy) , AF is the excess Helmholtz free 
energy (AF = AO, + p,(A(j>)/a), and AiV is the total number of molecules in a micelle. 

We briefly discuss two problems that can be encountered when calculating p m . First, all the densities are given 
in units of po/N, and all free energy densities in units of kTpo/N; on the other hand, a well know feature of the 
Gaussian model is that the spatial variation of the profiles is set by the combination a^/N. This means that, in order 
to obtain Afl/kT, we need the ratio between the two measures of volume, A = (a^/N) 3 /(N/ po) = poa 3 y^N', i.e., how 
many polymers would there be in a volume (a^/N) 3 . This is basically the parameter that describes the degree of 
concentration of a polymer solution, and by definition, A ^> 1 for a melt, but this value must be provided for each 
particular system. (Remarkably, our expression for v m does not depend on A.) 

Second, the calculation of AN is not obvious. There is always some arbitrariness in defining this magnitude; for 
copolymers, it is natural to propose AN C = A|A0|. But the micelles also include a certain number of solvent molecules 
(specially in the corona region). Looking at the profiles in Fig. 2, this would be the integral of the homopolymer 
profiles, dashed lines (together, of course, with a factor of 47rr 2 ), but restricted to the micellar region, since this 
integral diverges, unlike the one for copolymers. We propose to use a weight w(\<f>(r) — <fi b \), where w(x) is a function 
with a limit w(x) — > x for small x but w(x) — > 1 for large x to ensure the integration is properly restricted. 

Finally, we would like to comment that we have compared our calculation with a previous lattice model^. The 
authors use a form of SCFT, and Eq. J3J (including a proposal for v m ). The systems is slightly different, since they 
consider a simple solvent, not a homopolymer one, but the proper theory is very similar to the one written above: all 
the equations remain the same, except that the diffusion equations become simple Boltzmann distributions, since a 
solvent molecule has no internal conformations: 4>a,s(x) — oizq s (r) and q s (r) = exp[— awA(r)]. The appropriate choice 
of our parameters to mimic the lattice model is: / = 1/2, a = 1/100, \N — 124, and A = 10. Our results deviate 
markedly from this Reference: the location of the CMC is quite close, but in our case the volume fraction of micelles 
is high as soon as they become stable, so the assumption of isolated micelles breaks down, and one should consider 
instead the dense, periodic phase (very likely, a cubic phase) as the proper equilibrium structure. This difference 
is likely due to the lattice description of chain conformations in Ref. Il4l which can lead to an underestimation of 
the corresponding entropy and hence to a prediction of a much smaller micellar concentration. In addition, we have 
checked that spherical micelles are indeed the stable structure, by comparing with cylindrical ones and planar bilayers. 
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